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ABSTRACT 

Groundwater  storage  becomes  an  important  phenomenon  in  the  current  century  throughout  the  entire  globe. 
Basically,  groundwater  is  the  water  stored  deep  beneath  the  earth's  surface  in  underground  aquifers  and  it  is  less  vulnerable 
to  pollution  than  surface  water.  Thus,  it  is  an  essential  source  for  drinking  water.  But  today,  severe  decline  in  water  table 
due  to  anthropogenic  activities  becomes  a critical  problem  particularly  in  the  metropolitan  cities  with  high  population 
density.  This  water  storage  capacity  primarily  depends  on  the  balance  between  precipitaion  and  evapotranspiration  (ET). 
Various  meteorological  parameters  like  temperature,  wind  speed,  humidity,  solar  radiation  etc.  have  considerable  impact 
on  the  groundwater  storage.  It  is  extremely  necessary  to  understand  the  ET  of  the  region  before  town  and  country  planning. 
Rainwater  harvesting  is  a common  practice  for  groundwater  recharge.  In  urban  areas,  water  bodies  like  ponds  and  lakes 
may  play  a crucial  role  for  groundwater  storage.  The  present  study  focused  on  the  role  played  by  a pond  as  natural  water 
harvesting  system  in  an  academic  institution  in  Kolkata.  EA056  Penman-Monteith  method  has  been  extensively  used  for 
estimating  reference  evapotranspiration  for  the  study  area  from  201  Ito  2016.  The  results  indicate  higher  ETq  during  the 
summer  season  while  higher  water  storage  in  the  monsoon  season.  The  months  from  June  to  September  have  considerably 
high  water  storage.  47.16  % of  the  total  annual  rainfall  has  been  stored  in  the  pond  for  the  last  five  years  continuously.  It 
reflects  a great  potential  of  groundwater  storage  using  natural  rainwater  harvesting  system  in  the  adjacent  locality. 

KEYWORDS:  EA056  Penman-Monteith  Method,  Groundwater  Storage,  Meteorological  Parameters,  Rainwater 
Harvesting,  Reference  Evapotranspiration 

INTRODUCTION 

Groundwater  is  recharged  from,  and  eventually  flows  to,  the  surface  naturally;  natural  discharge  often  occurs 
at  springs  and  seeps,  and  can  form  oases  or  wetlands.  Groundwater  is  also  often  withdrawn  for  agricultural,  municipal, 
and  industrial  use  by  constructing  and  operating  extraction  wells.  Groundwater  primarily  depends  on  precipitation  and 
evapotranspiration.  It  is  directly  related  to  the  occurrence  of  rainfall  and  inversely  related  to  the  evaporation  and  surface 
runoff  Geology,  soil  and  natural  vegetation  play  an  important  role  in  groundwater  storage.  Permeable  rock  with  water 
retention  capacity  or  aquifer  is  essential  for  groundwater  recharge.  Clay  soil  with  high  porosity  provides  obstruction  in 
groundwater  recharge.  Grassland  region  has  more  potential  for  groundwater  storage  than  dense  forest  due  to  its  low  rate  of 
evapotranspiration. 

Evaporation  comprises  an  important  part  of  the  thermal  balance  occurring  on  the  earth  surface.  It  is  also  a part  of 
water  budget  and  the  surface  heat  and  water  conditions  are  the  determinants  of  formation  of  the  surface  ecological 
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environment.  Thus  the  study  of  evaporative  capacity  of  the  land  is  always  one  of  the  major  problems  in  hydrology  and 
geosciences  (Ping  et  ah,  2009).  The  rate  of  reference  evapotranspiration  is  generally  calculated  from  a reference  surface 
where  it  is  denoted  as  ETq  but  adequate  amount  of  water  is  required  there.  Reference  surface  for  the  purpose  is  a 
hypothetical  grass  which  is  a reference  crop  with  specific  characteristics.  With  the  help  of  ETq,  many  research  have  been 
done  like  aridity/humidity  conditions  (Wu  et  ah,  2006),  ecosystem  models  (Eisher  et  ah,  2005),  estimation  of 
rainfall-runoff  and  water  use  in  agriculture  (Allen,  2000;  Hunsaker  et  ah,  2002).  The  evapotranspiration  of  reference  crop 
is  found  on  the  basis  of  the  meteorological  data  and  the  calculation  is  given  by  EAO  depending  on  the  meteorological 
parameters  and  meteorological  data. 

Thus  it  is  difficult  to  apply  the  method  in  the  areas  where  there  is  lack  of  meteorological  observation  data,  and 
recommended  the  use  of  evaporating  pan  observation  data  for  determining  the  reference  crop  evapotranspiration  by  EAO 
(Allen  et  ah,  1998).  Various  methods  for  empirical  estimation  of  ETq  requires  proper  and  accurate  measurements  of 
different  parameters  like  temperature,  humidity,  wind  speed,  sunshine,  solar  radiation  etc.  But  measurement  of  all  these 
parameters  in  a place  is  quite  few  particularly  in  any  developing  country.  There  are  also  local  changes  in  ETq  on  the  basis 
of  distance  from  the  weather  station  (Hubbard,  1994;  Pielke  et  al.,  2000)  and  proper  integration  of  different  climatic 
parameters  are  required  which  also  affect  the  accuracy  of  ETq  (Meek  and  Hatfield,  1994;  Allen,  1996).  Various  methods 
for  calculations  are  there  among  which  pan  evaporimeter  and  some  ET  models  use  only  temperature  and  thus  are  less 
complex  (Magliulo  et  al.,  2003).  The  Hargreaves  equation  also  needs  only  minimum  and  maximum  temperature 
(Hargreaves  and  Samani,  1985)  and  extraterrestrial  radiation  (Droogers  and  Allen,  2002).  Thus  ETq  can  be  measured  by 
various  models  and  methods  by  different  weather  parameters  (Thornthwaite,  1948;  Penman,  1948;  Priestley  and  Taylor, 
1972;  Hargreaves,  1994;  Hargreaves  and  Allen,  2003).  Allen  et  al.,  (1998)  developed  a model  for  ETq  which  was  published 
by  the  Eood  and  Agriculture  Organization  of  the  United  Nations  as  Penman-Monteith  (EA056-PM)  by  using  a hypothetic 
reference  crop  which  is  approved  for  the  arid  and  humid  climatic  conditions.  A number  of  research  works  and  applications 
on  natural  rainwater  harvesting  have  been  experienced  in  parts  of  South  Indian  states.  A pond  and  its  surroundings  were 
considered  for  the  cutTent  work  which  may  be  treated  as  a model  study  for  further  application.  Present  study  involves  the 
groundwater  storage  estimation  through  ETq  computation  and  analysis  of  variation  for  the  combined  area  of  pond,  play 
ground  and  shaded  portion  in  the  Dinabandhu  Andrews  College  premises  with  parameters  like  minimum  and  maximum 
temperature,  humidity,  wind  speed,  sunshine  hours  and  solar  radiation  for  the  last  five  years  (201 1 to  2016). 

STUDY  AREA 

The  study  area  is  a pond  lies  in  Dinabandhu  Andrews  College  premises,  Garia,  Kolkata,  India  which  extends  from 
22°28'04"N  to  22°28’06"N  latitude  and  from  88°22'41"N  to  88“22'43"E  longitude  with  an  area  of  1459.92  sq  m.  The  total 
area  including  the  shaded  area  and  the  ground  area  used  to  estimate  the  recharge  the  groundwater  level  beneath  the  surface 
of  the  pond  is  4600  sq  m.  The  climate  is  mostly  humid  and  witnesses  an  annual  average  rainfall  of  140.27  mm.  The 
location  map  of  the  study  area  is  presented  in  Eigure  1 . 

A number  of  satellite  data  acquired  to  show  the  seasonal  synoptic  views  of  the  pond  and  the  adjacent  parts  of  the 
Dinabandhu  Andrews  College  Premises  (Eigure  2).  The  shaded  area  is  made  of  asbestos  and  it  is  slightly  inclined  to  the 
pond.  Thus,  there  is  an  artificial  runoff  process  observed  over  the  asbestos  shade.  The  ground  is  partially  covered  with 
grass  and  bordered  by  a number  of  trees.  Moreover,  the  ground  has  a very  gentle  slope  to  the  pond  and  there  is  a channels 
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between  the  ground  and  the  pond.  Hence,  the  surface  runoff  water  flows  over  the  play  ground  is  normally  accumulated  into 
the  pond.  Besides,  the  pond  has  no  significant  outlet  or  no  domestic  use. 


Figurel:  Location  of  the  Study  Area 


Figure  2:  Multi-Temporal  Satellite  Imageries  of  Different  Seasons  Showing 
Pond,  Play  Ground  and  Shaded  Area  in  Dinahandhu  Andrews  College  Premises 

MATERIALS  AND  METHODS 

Estimation  of  ETo  with  Penman-Monteith  Equation 

For  calculating  the  ETq,  climatic  station  in  Kolkata  was  used  in  the  study.  Different  climatic  variables  like  rainfall, 
minimum  and  maximum  air  temperature,  relative  humidity,  air  pressure,  wind  speed  and  solar  radiation  of  last  five  years 
(2011  to  2016)  were  considered.  Average  monthly  datasets  have  been  used  for  estimating  the  average  monthly  ETO  taken 
from  the  Water  Development  and  Management  Unit  and  the  Climate  Change  and  Bio-energy  Unit  of  EAO,  2006. 

The  Penman-Monteith  model  was  formulated  on  the  basis  of  hypothetical  green  grass  reference  surface  in  which 
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case  the  height  of  the  grass  is  presumed  to  be  0.12m  with  a surface  resistance  of  70s  m-1  and  with  an  albedo  of  0.23 
(Allen  et  al.,  1998).  It  has  been  approved  by  Food  and  Agriculture  Organization  (FAO-  56).  Penman  Monteith  equations 
are  given  below  in  sequence: 


0.40^(r,  - g)+  (e,  - e J 


ET„=- 


T + 273 


A+y(l+0.34U2) 


(1) 


Where, 

ETo  = reference  evapotranspiration  [mm  day'*] 

R„  = net  radiation  at  the  crop  surface  [MJ  m'^  day'*] 

G = soil  heat  flux  density  [MJ  m'^  day'*] 

T = mean  daily  air  temperature  at  2 m height  [°C] 

U2  = wind  speed  at  2 m height  [m  s'*] 
es  = saturation  vapour  pressure  [kPa] 

ea  = actual  vapour  pressure  [kPa]  es  - e^  = saturation  vapour  pressure  deficit  [kPa] 

A = slope  vapour  pressure  curve  [kPa  °C'*] 
y = psychrometric  constant  [kPa  °C'*] 

ETq=  reference  evapotranspiration  (mmd-1) 

Slope  of  Saturation  Vapour  Pressure  (A) 

The  slope  of  the  relation  between  saturation  vapour  pressure  and  temperature  A is  required  for  calculating 
reference  evapotranspiration.  At  a given  temperature,  the  slope  is  given  by: 


2504expr^Z^l 

4098^(t)  _ l^T+  237.2j 

(T+  237.3)^  (T+  237.3)^ 

Where, 

A = slope  vapour  pressure  curve  [kPa  C'*] 

T = air  temperature  [°C] 

e°  (T)  = saturation  vapour  pressure  at  temperature  T [kPa] 

Atmospheric  Pressure  (P) 


(2) 


(3) 


Where, 
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P = atmospheric  pressure  (KPa) 


Z = elevation  above  sea  level  (m) 


Psychrometric  Constant  (y) 


(4) 


Where, 

y=  psychrometric  constant  [kPa  C'*] 

Cp  = specific  heat  of  moist  air  = 1.013  [kJ  kg'*  °C'*] 

P = atmospheric  pressure  [kPa] 

e=  ratio  molecular  weight  of  water  vapour/dry  air  = 0.622 
h=  latent  heat  of  vaporization  [MJ  kg'*] 

Saturation  Vapour  Pressure  (es) 

As  saturation  vapour  pressure  is  related  to  air  temperature,  it  can  be  calculated  from  the  air  temperature.  The 
relationship  is  expressed  by: 


Where, 

e°(T)  = saturation  vapour  pressure  function  [kPa] 

T = air  temperature  [°C] 

Due  to  the  non-linearity  of  the  above  equation,  the  mean  saturation  vapour  pressure  for  a day,  week,  decade  or 
month  should  be  computed  as  the  mean  between  the  saturation  vapour  pressure  at  the  mean  daily  maximum  and  minimum 
air  temperatures  for  that  period: 


(5) 


e°(T,J+e°(T,J 


2 


(6) 


3.1.5  Actual  Vapour  Pressure  (ea) 


2 


(7) 


Where, 


ea  = actual  vapour  pressure  [kPa] 


e°(T„,in)  = saturation  vapour  pressure  at  daily  minimum  temperature  [kPa] 
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e°(Tniax)  = saturation  vapour  pressure  at  daily  maximum  temperature  [kPa] 

RHmax  = maximum  relative  humidity  [%] 

= minimum  relative  humidity  [%] 

Extraterrestrial  Radiation  for  Daily  Periods  (Ra) 

The  extraterrestrial  radiation,  R^,  for  each  day  of  the  year  and  for  different  latitudes  can  be  estimated  from  the 
solar  constant,  the  solar  declination  and  the  time  of  the  year  by; 


Where, 

Ra  = extraterrestrial  radiation  [MJ  m'^  day'*] 

Gsc  = solar  constant  = 0.0820  MJ  m'^  min'* 
dr  = inverse  relative  distance  Earth-Sun 
w s = sunset  hour  angle  [rad] 
cp  = latitude  [rad] 

6=  solar  decimation  [rad] 

Ra  is  expressed  in  the  above  equation  in  MJ  m'^  day  *.  The  corresponding  equivalent  evaporation  in  mm  day  * is 
obtained  by  multiplying  Ra  by  0.408.  The  latitude,  cp,  expressed  in  radians  is  positive  for  the  northern  hemisphere  and 
negative  for  the  southern  hemisphere. 

Solar  Radiation  (Rj) 

If  the  solar  radiation,  R,,  is  not  measured,  it  can  be  calculated  with  the  Angstrom  formula  which  relates  solar 
radiation  to  extraterrestrial  radiation  and  relative  sunshine  duration: 


Where, 

R,  = solar  or  shortwave  radiation  [MJ  m'^  day'*] 
n = actual  duration  of  sunshine  [hour] 

N=  maximum  possible  duration  of  sunshine  or  daylight  hours  [hour] 
n/N  = relative  sunshine  duration  [-] 

Ra  = extraterrestrial  radiation  [MJ  m'^  day'*] 

as  = regression  constant,  expressing  the  fraction  of  extraterrestrial  radiation  reaching  the  earth  on  overcast  days 


(8) 


(9) 
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(n  = 0) 

as+bs  = fraction  of  extraterrestrial  radiation  reaching  the  earth  on  clear  days  (n  = N) 

Clear-Sky  Solar  Radiation  (Rs„) 

The  calculation  of  the  clear-sky  radiation,  Rso,  when  n = N,  is  required  for  computing  net  longwave  radiation. 

Rso  = (as+bJRa  (10) 

When  calibrated  values  for  Uj  and  b^  are  available 

R,o  = (0.75  + 2 lO'^z)  Ra  (11) 

When  calibrated  values  for  Uj  and  b^  are  not  available 

3.1.9  Net  Solar  or  Net  Shortwave  Radiation  (R„s) 

The  net  shortwave  radiation  resulting  from  the  balance  between  incoming  and  reflected  solar  radiation  is  given 
by: 

R„,  = (l-a)Rs  (12) 

Where, 

R„s  = net  solar  or  shortwave  radiation  [MJ  m'^  day'*] 

a = albedo  or  canopy  reflection  coefficient,  which  is  0.23  for  the  hypothetical  grass  reference  crop 
[dimensionless] 

R,  = the  incoming  solar  radiation  [MJ  m'^  day'*] 

Net  Longwave  Radiation  (R„i) 

The  rate  of  longwave  energy  emission  is  proportional  to  the  absolute  temperature  of  the  surface  raised  to  the 
fourth  power.  This  relation  is  expressed  quantitatively  by  the  Stefan-Boltzmann  law.  The  net  energy  flux  leaving  the  earth's 
surface  is,  however,  less  than  that  emitted  and  given  by  the  Stefan-Boltzmann  law  due  to  the  absorption  and  downward 
radiation  from  the  sky.  Water  vapour,  clouds,  carbon  dioxide  and  dust  are  absorbers  and  emitters  of  longwave  radiation. 
Their  concentrations  should  be  known  when  assessing  the  net  outgoing  flux.  As  humidity  and  cloudiness  play  an  important 
role,  the  Stefan-Boltzmann  law  is  corrected  by  these  two  factors  when  estimating  - the  net  outgoing  flux  of  longwave 
radiation.  It  is  thereby  assumed  that  the  concentrations  of  the  other  absorbers  are  constant: 


"*"  ^mintd 


[0.34-0.14ji7ll.35^-0.35 


Where, 

Rni  = net  outgoing  longwave  radiation  [MJ  m'^  day'*] 
s = Stefan-Boltzmann  constant  [4.903  10'^  MJ  K'"*  m'^  day'*] 


(13) 


Impact  Factor(JCC):  2.7341  - This  article  can  be  downloaded  from  www.impactiournals.us 


106 


Somnath  Mukhopadhyay  & Subhanil  Guha 


Tmax,  K = maximum  absolute  temperature  during  the  24-hour  period  [K  = °C  -H  273. 16] 

Tmin,  K = minimum  absolute  temperature  during  the  24-hour  period  [K  = °C  H-  273.16] 
e^  = actual  vapour  pressure  [kPa] 

Rs/Rso  = relative  shortwave  radiation  (limited  to  < 1.0) 

Rs  measured  or  calculated.  (Equation  35)  solar  radiation  [MJ  m'^  day'*] 

Rso  calculated  (Equation  36  or  37)  clear-sky  radiation  [MJ  m'^  day'*] 

Net  Radiation  (Rn) 

The  net  radiation  (R„)  is  the  difference  between  the  incoming  net  shortwave  radiation  (R„s)  and  the  outgoing  net 
longwave  radiation  (R„i): 

Rn  = Rn.s  - Rnl  (14) 

Where, 

Rns  = incoming  net  shortwave  radiation 
Rnl  = outgoing  net  longwave  radiation. 

Wind  Speed  (U2) 

Eor  the  adjustment  of  the  wind  speed  data  obtained  from  instrument  which  is  placed  at  elevations  other  than  the 
standard  height  of  2 m,  following  calculation  is  done: 


4.87 

ln(67.8z-5.42) 

Where, 

U2  = wind  speed  at  2 m above  ground  surface  [m  s'*], 

uz  = measured  wind  speed  at  z m above  ground  surface  [m  s'*], 

z = height  of  measurement  above  ground  surface  [m]. 

Soil  Heat  Flux  (G) 

Complex  models  are  available  to  describe  soil  heat  flux.  Because  soil  heat  flux  is  small  compared  to  R„, 
particularly  when  the  surface  is  covered  by  vegetation  and  calculation  time  steps  are  24  hours  or  longer,  a simple 
calculation  procedure  is  presented  here  for  long  time  steps,  based  on  the  idea  that  the  soil  temperature  follows  air 
temperature: 


G = c, 

' At 


(16) 


Where, 
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G = soil  heat  flux  [MJ  m'^  day'*] 

Cs  = soil  heat  capacity  [MJ  m'^  °C'*] 

Ti  = air  temperature  at  time  i [°C] 

Ti_i  = air  temperature  at  time  i-1  [°C] 

A t = length  of  time  interval  [day] 

A z = effective  soil  depth  [m] 

As  the  soil  temperature  lags  air  temperature,  the  average  temperature  for  a period  should  be  considered  when 
assessing  me  daily  soil  heat  flux,  i.e.,  A t should  exceed  one  day.  The  depth  of  penetration  of  the  temperature  wave  is 
determined  by  the  length  of  the  time  interval.  The  effective  soil  depth,  A z,  is  only  0.10-0.20  m for  a time  interval  of  one  or 
a few  days  but  might  be  2 m or  more  for  monthly  periods.  The  soil  heat  capacity  is  related  to  its  mineral  composition  and 
water  content. 

RESULTS  & DISCUSSIONS 

Different  climatic  variables  of  minimum  and  maximum  temperature,  humidity,  wind  speed,  sunshine  duration  and 
radiation  have  been  shown  for  5 years  from  2011  to  2016.  The  minimum  temperature  is  highest  in  May-June  up  to 
September  after  which  it  started  to  decrease.  The  minimum  temperature  varies  from  26.5°C  in  May)  to  13.8°C  in  January. 
The  average  monthly  maximum  temperature  is  highest  in  the  month  of  April  and  May  with  more  than  35**C  and  lowest  in 
January  and  December  (<  27°C).  The  relative  humidity  is  maximum  from  June  to  September  (>  90%)  while  lowest 
humidity  is  found  from  February  to  March  (<  66%).  Wind  speed  is  maximum  from  April  to  July.  It  varies  from  69.1 
km/day  in  October- January  to  190.1km/day  during  April-May.  The  maximum  and  minimum  sunshine  hours  varies  from 
8.19  hr/day  in  April  when  sky  remains  very  clear  to  less  than  4.8  hr/day  during  the  monsoon  period.  The  solar  radiation  is 
highest  in  April  of  around  21.8  MJm'^d  * to  around  14.85  MJm'M  * during  December.  Rainfall  varies  from  5.32  mm  in 
December  to  357.55  mm  in  August.  Reference  evapotranspiration  (ETq)  has  been  calculated  with  the  Penman-Monteith 
method  in  mm/day.  The  value  of  ETq  is  ranging  between  2.52  in  December  to  5.24  in  April  and  May  (Eigure  3). 

The  mean  monthly  ETq  was  calculated  with  help  of  EAO  56  Penman-Monteith  method  to  show  the  variation  of 
ETo  in  the  study  area  which  is  an  extremely  important  variable  to  be  considered  for  groundwater  storage.  The  result 
indicates  that  ETq  is  highest  in  the  month  of  May  (5.24  mm  day'*)  to  lowest  in  the  month  of  December  (2.52  mm  day'*). 
Mean  monthly  ETq  of  the  area  is  3.68  mm  day  *.  To  calculate  the  total  water  loss  due  to  evapotranspiration,  impact  of  crop 
coefficient  (Kc)  should  be  considered  which  varies  with  different  land  use  and  land  cover  types. 

In  the  crop  coefficient  approach  the  crop  evapotranspiration  (ET^)  is  calculated  by  multiplying  the  reference  crop 
evapotranspiration,  (ETo)  by  a crop  coefficient  (K^): 

ETo  = Ko  X ETo  (17) 

Where, 

ETo  = crop  evapotranspiration  [mm  d'*] 
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K(.  = crop  coefficient  [dimensionless] 

ETo  = reference  crop  evapotranspiration  [mm  d'*] 

The  soil  surface  wetness  and  the  fraction  of  ground  covered  by  vegetation  influence  the  surface  resistance,  r^. 
Following  soil  wetting,  the  vapour  transfer  rate  from  the  soil  is  high,  especially  for  crops  having  incomplete  ground  cover. 
The  combined  surface  resistance  of  the  canopy  and  of  the  soil  determines  the  (bulk)  surface  resistance,  rs.  The  surface 
resistance  term  in  the  Penman-Monteith  equation  represents  the  resistance  to  vapour  flow  from  within  plant  leaves  and 
from  beneath  the  soil  surface. 

In  the  present  research  work,  values  of  K;.  are  considered  as  1.05  for  open  water  source  of  the  pond,  0.4  for  silty 
clay  soil  surface  of  the  play  ground  and  O.lfor  the  shaded  area  in  the  college  premises.  The  annual  actual  water  loss 
through  evapotranspiration  is  estimated  as  3618.72  cubic  metre[  (2056.9  cubic  metre  from  pond)  + (1517.38  cubic  metre 
from  play  ground)  + (44.44  cubic  metre  from  shaded  area)]. 

The  play  ground  has  mainly  silty  clay  type  of  soil  texture.  It  shows  a low  rate  of  infiltration  (5-10  mm/hour) 
during  heavy  precipitation.  This  type  of  soil  allows  high  rate  of  surface  runoff.  Apart  from  that,  the  ground  has  numerous 
channels  through  which  water  can  move  into  the  pond.  Actually,  this  ground  is  considered  as  an  inlet  of  the  pond. 
Approximate  10%  of  the  total  volume  of  rainfall  is  infiltrated  through  the  play  ground. 

The  total  volume  of  groundwater  storage  in  the  pond  is  computed  as  3651.37  cubic  metre  per  year.  It  means  that 
the  study  area  is  annually  recharged  with  47.16%  of  the  total  annual  precipitation.  The  results  vividly  indicate  that  the  pond 
of  the  college  premises  may  be  considered  as  a highly  effective  groundwater  recharge  potential  zone. 


Figure  3:  Pond  Acts  as  a Natural  Water  Harvesting  System  (Dinabandhu  Andrews  College,  Kolkata) 
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Figure  4:  Average  Monthly  Climatic  Variables  (2011-2016) 


CONCLUSIONS 


The  study  involves  the  computation  and  analysis  of  one  of  the  most  important  climatic  variable,  i.e. 
evapotranspiration  for  the  computation  of  annual  groundwater  storage.  This  study  shows  the  monthly  ET  variation  which  is 
very  useful  for  analysis  of  variation  in  land  use  and  land  cover  planning.  The  disparity  of  ET  in  different  parts  of  the  area 
indicates  the  potential  of  water  storage  by  means  of  proper  planning  and  conservation  of  water  resources.  It  can  solve  the 
problem  of  water  scarcity  of  the  locality  in  the  very  near  future.  Dinabandhu  Andrews  College  plays  a major  role  in  the 
groundwater  storage  through  natural  water  harvesting  technique  and  its  proper  conservation.  EAO-56  Penman-Monteith 
method  is  a globally  accepted  and  extremely  important  empirical  model  to  calculate  ETq  which  is  also  an  important 
meteorological  parameter.  This  research  may  be  considered  as  a positive  outset  for  the  entire  locality  and  the  city  to 
groundwater  storage  through  natural  harvesting. 
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